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Quantum optics with 
nitrogen-vacancy centers in diamond 


1.1 Background 


The study of classical mechanics, thermodynamics, and electromagnetism have brought 
on great technological advances over the course of history. We now live in an exciting 
era where quantum mechanics begins to join these other fields of physics in enabling 
new technologies with unprecedented capabilities. Theoretical and experimental ef¬ 
forts have moved toward harnessing the power of various different physical systems for 
applications in metrology, quantum information, and furthering our understanding of 
quantum mechanics itself. 

Some of the most successful of applications of quantum mechanical systems to date 
have been in the field of metrology, or the accurate sensing and measurement of forces. 


fields, frequencies, etc. (Chou et al.^ 2010 Nicholson et al.^ 2012 Benko et a/., 2012 


Kessler et a/., 2012). At the same time, the field of quantum information and commu¬ 


nications has become the focus of much effort in recent years. There has been a large 
number of theoretical proposals for using quantum systems for efficient computation. 


simulation of complex systems, secure dissemination of information, etc. (Shor, 1994 


Kimble, 2008 Lloyd, 1996). However, the practical implementation of these propos¬ 


als in real physical systems remains an open challenge. A variety of different avenues 
are being actively explored, ranging from photons, neutral atoms, ions, impurities 
and structures in solid state systems, to more exotic materials such as electrons on 


the surface of liquid helium (Monroe and Kim, 2013 Awschalom et a/., 2013 Clarke 


and Wilhelm, 2008 Stern and Lindner, 2013 Platzman and Dykman, 1999). A com¬ 


mon theme for the implementation of quantum technologies involves addressing the 
seemingly contradictory need for controllability and isolation from external effects. 
Undesirable effects of the environment must be minimized, while at the same time 
techniques and tools must be developed that enable us to interact with the system in 
a controllable and well-defined manner. These lectures address several aspects of this 
theme with regards to a particularly promising candidate for developing applications 
in both metrology and quantum information: the nitrogen-vacancy (NV) center in di¬ 
amond. In particular, we explore how the NV center can be investigated, controlled, 
and efficiently coupled with optical photons, and tackle the challenges of controlling 
the optical properties of these emitters inside a complex solid-state environment. Some 


recent reviews on NV centers can be found in (Jelezko and Wrachtrup, 2006; Doherty 


et al.^ 2013 Doherty et al.^ 2012). 





































2 Quantum optics with nitrogen-vacancy centers in diamond 


As with other solid state systems such as quantum dots, rare earth dopants in 
inorganic crystals, and superconducting qubits, NV centers present the attractive fea¬ 
ture of eliminating the need for trapping and cooling, as is usually required for atomic 
systems. In addition to greatly reducing the complexity of experiments, solid state 
systems offer the possibility of straightforward integration with fabricated structures 
and development of scalable quantum devices. The issue of scalability, however, is not 
simply a matter of being able to create a large number of physical devices in a small 
volume. In order to truly move beyond proof-of-principle experiments, the properties 
of many qubits and devices must all be well controlled in a reproducible, reliable, and 
efficient manner. In other words, while it is important to demonstrate that a certain 
quantum operation, coherence time, metrological sensitivity, etc. is possible in one 
particular device with favorable characteristics, it must also be possible to reproduce 
these characteristics for a large number of systems. 

Diamond is a particularly promising host material for optically active quantum 
emitters due to its wide bandgap (5.5 eV) and large transparency window. It supports 
a huge variety of impurities and defects that give rise to optical bands at energies 
throughout the infrared and visible spectrum (Zaitsev, 2010). The NV center is unique 
in that it combines well-defined optical transitions with long-lived spin sub-levels in 
the ground state, much like alkali atoms with transitions involving hyperfine sublevels. 
The present lectures focus on quantum optical experiments with NV centers. In recent 
years, a comprehensive theoretical model has been developed to describe the spin and 


optical properties of the NV center, which has been reviewed in (Maze et a/., 2011 


Doherty et a/., 2012 Doherty et a/., 2010). We begin by summarizing this theoretical 


framework, which allows us to understand the underlying physics of the subsequent 
experimental results. We then describe the experimental techniques that are used to in¬ 
vestigate and address the optical properties of NV centers. As part of this description, 
we will encounter some general characteristics of the system that will be come impor¬ 
tant issues to be addressed later, such as the Debye-Waller factor and the sensitivity 
of the excited state to the local environment. Based on our theoretical understanding 
and spectroscopic results, we show that NV centers offer a rich set of possibilities for 
experiments demonstrating coherent optical control of the system, several of which are 
summarized in the remainder of the chapter. 


1.2 Level structure and polarization properties of the NV center 


As is the case in any system, the electronic structure of the NV center can be analyzed 
by considering its symmetry, in addition to the various interactions that are present. 
Our description below starts with electronic states that obey the Csy symmetry of 
the NV center. We then consider the effect of various interactions one by one, starting 
from the the Coulomb interaction, which sets the energy scale of the optical transitions. 
Spin-orbit and spin-spin interactions act to lift the degeneracy of the different states 
in the ground and excited state manifolds. Finally, we consider the effect of external 
perturbations such as electric and magnetic fields that change the nature and energies 
of the unperturbed states. More detailed analysis of many of the concepts presented 
in this section can be found in (Maze et a/., 2011). 
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Fig. 1.1 Single electron orbitals of the NV center. The nitrogen and three carbons surround¬ 
ing the vacancy are shown, and the perspective is along the NV axis. The color scale roughly 
represents the sign and occupation of each orbital. 

1.2.1 Electronic states of the NV center 


The negatively-charged NV center has six electrons, five of which are from the nitro¬ 
gen and the three carbons surrounding the vacancy. They occupy the orbitals states 
a'l, ai, e^, which are combinations of the four dangling bonds surrounding the va¬ 
cancy that satisfy the symmetry imposed by the nuclear potential. They can be written 
as 


a'l = a 
ai = p 
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where ( 7 i, ( 72 , (73 and (Jn are the dangling bonds from the three carbons and the 
nitrogen, respectively. These orbitals transform as the irreducible representations of 


the group and have been extensively used by many authors (Loubser and Wyk, 


1978 Goss et a/., 1996 

1 Lenef and Rand, 1996 Gali et al.^ 2008). They are schematically 

illustrated in Figure | 

1 . 1 | 

to give an idea of their symmetries. From these pictorial 


representations, one can deduce that, for example, {ey\y • r\ex) = {ex\x • r\ex) = 
— {ey\x ' r\ey) 7 ^ 0. This means that the je^,) and |e^) states have permanent electric 
dipole moments, which will become important later for understanding the electric 
field sensitivity of the NV center’s excited states and spectral diffusion of the optical 
transitions. In addition, one also sees that {a\x • r\ex) = (a| ^ • r je^) 7 ^ 0, which gives 
rise to the transition dipole moments between the ground and excited states. These 
properties will be discussed in more detail in the following sections. 

In the ground state of the NV center, the a'^ and ai states, which are lowest 
in energy, are filled by four electrons, as shown in Figure | 1 . 2 | The remaining two 
electrons occupy the degenerate orbitals Cx and Cy. The orbitals Cx and Cy can be 
viewed as p-type orbitals and e+ = — e_ = Cx — iCy are analogous to p 

states with definite orbital angular momentum. An antisymmetric combination of the 
orbital states minimizes the Coulomb energy and results in a spin-triplet ground state 
manifold 
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Fig. 1.2 Ground and excited states associated with the NV center’s optical transitions. The 
four symmetrized states are filled by six electrons, which can also be viewed as two holes 
(indicated with green arrows). The Coulomb interaction defines the optical transition energy 
of the NV center. The effect of the spin-orbit and spin-spin interactions are indicated for both 
the ground and excited states. In the ground state, the application of a magnetic field further 
splits the \ms = if) states due to the Zeeman interaction. In the excited state, crystal strain 
modifies the energies of sublevels. See (Maze et ah, 2011) for details. 
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where |±1) , |0) correspond to the rris = ±1,0 states, respectively. The notation |^o) = 
— CyCx) indicates that the orbital state has 0 orbital angular momentum pro¬ 
jection along the NV axis. A 2 denotes the orbital symmetry of the state, which is 
determined by the symmetries of the and Cy orbitals. From here on, we will often 
denote these ground states simply as |0) and |±1) when it’s clear that we are referring 
to the full electronic state and not just the spin projections. 

The relevant excited state for the optical transitions of the NV center is a pair of 
triplets which arises from the promotion of one of the electrons occupying the orbital 
ai to the Cx or Cy orbitals (Goss et al.^ 1996). Note that this state can be modeled 


by one hole in the orbital e and another hole in the orbital ai, i.e. a triplet in the ae 
electronic configuration (similarly, the ground state can be modeled by two holes in the 
electronic configuration). A total of six states can be formed in this configuration 


and their symmetries are determined by a group theoretical analysis (Maze et al.^ 


2011 ) 
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1 ^ 2 ) = |-E_) 0 1+1) + \E+) 0 1-1) 

\E,) = \X) 0 |0) 

\Ey) = \Y) 0 |0) (1.6) 

|Ei) = |.B_) 0 1-1) - |E+) 0 1+1) 

\E2) = \E_) 0 1-1) + \E+) 0 1+1) 


where we have named the first four states as Ai, A 2 , and Ey according to their 
symmetries and named the last two states as Ei and E 2 since they also transform 
according to the irreducible representation E. Here, \E±) = |ae± — e±a) and \X{Y)) = 

As spin-orbit and spin-spin interactions are invariant under any operation of the 
Csy group, the states given in equation (1.6) are eigenstates of the full Hamiltonian 
including these interactions and in the absence of any perturbation such as magnetic 
field and/or crystal strain. The Hamiltonians for the spin-orbit and spin-spin interac¬ 
tions can be written in the basis of the six states listed above as 


Hso = Xz + |^2)(A2| - - |^2)(^2|) (1.7) 

Hss = A (|Ai)(Ai| + |A2)(A2| + + \E2){E2\) 

-2A {\E,){E,\ + \Ey){Ey\) + A' {\A2){A2\ - |^i)(^i|) 

+A" {\E^){Ey\ + \Ey){E^\- i\E2){E,\ + i\E,){E2\) (1.8) 


where A;^ is the axial spin-orbit interaction and 3A 1.42 GHz and A' « 1.55 GHz 
characterize the spin-spin induced zero-field splittings ( |Lenef and Rand, 1996[ [Rogers 


et a/., 2009 Maze et a/., 2011). It can be seen that spin-orbit interaction splits states 


with different total angular momentum, i.e., the pairs (^ 1 ,^ 2 ), (E^^Ey) and (^ 1 ,^ 2 ) 


split from each other by about 5.5 GHz ( Lenef and Rand, 1996| Manson et a/., 200^ 
Batalov et a/., 2009). The A" term leads to non-spin preserving cross transitions, which 


results to optical pumping of the spin states and plays a role in, for example, limiting 
the efficiency of spin readout using resonant excitation. The spin-spin interaction also 
plays a crucial role in the stability of the states \Ai) and IA 2 ). Thi s interaction shifts 
up the non-zero spin states (Ai, A 2 , Ei and E 2 ) by ^ GHz (Fuchs et al.^ 2008 


m c 

Pi, 


Batalov et a/., 2009| and shift the states \Ex) and 
the gap between 


\Ey) down by ^ |l.42 GHz. Thus, 


A 2 ) and (Ex^Ey) is increased, but the gap between the states 
{Ex, Ey) and (F^i, £^ 2 ) is reduced. In addition, the spin-spin interaction splits the states 
IA 2 ) and |Ai) by ^ 3.3 GHz ( [Batalov et al, 2009 ). Thus, for relatively low strain, we 
obtain a very robust IA 2 ) state with stable symmetry properties that are protected by 
an energy gap arising from the spin-orbit and spin-spin interactions. 

We note here that the zero-field splitting in the ground state also arises from the 
spin-spin interaction, and gives an energy difference of Ags ^2.88 GHz between the 
\ms = 0) and = ±1) states. 


1.2.2 Properties of optical transitions 

Once the wavefunctions are known, it is possible to calculate the selection rules of 
optical transitions between the triplet excited state and the triplet ground state. The 
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Table 1.1 Selection rules for optical transitions between the triplet excited state (ae) and 
the triplet ground state (e^). Linear polarizations are represented by x and y, while circular 
polarization are represented by (j± = x zL iy. As an example, a photon with cr+ polarization 
is emitted when the electron decays from state A 2 to state ^^ 2 -. 
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dipole moment between the ground and excited state is produced by the hole left 
in the a orbital under optical excitation. As mentioned earlier, the matrix elements 
{a\x ' f\ex) and {a\y • r\ey) are non-zero, where x and y represent the polarization of 
the involved photon. We can then calculate the selection rules for transitions between 
every pair of ground and excited states, as shown in Table [m As expected, these 
selection rules conserve the total angular momentum of the photon-NV center system. 

The allowed optical transitions and their polarization properties indicate several 
possible A schemes in the NV center. However, while these properties are relatively 
robust in the \Ai) and IA 2 ) states, the \Ai) state is coupled non-radiatively to a 
metastable singlet state, which then decays to the ground state |0). This results in 
leakage out of the A system consisting of the \Ai) state and the |±1) states. Thus, 
we find that the IA 2 ) state provides the an almost ideal closed A scheme for experi¬ 
ments such as spin-photon entanglement, as described in Section [T^ and ( |Togan et a/., 
2010). On the other hand, the open A system involving \ Ai) can be useful for efficient 
detection of dark states during coherent population trapping, as shown in Section [T77| 


and (Togan et a/., 2011) 


1.2.3 The effect of strain 

We now discuss the effect of local strain on the properties of the optical transitions in 
order to understand variations between different NV centers and deviations from the 
unperturbed system described above. This perturbation splits the degeneracy between 
the Cx and Cy orbitals and results in their mixing. In the limit of high strain (larger 
than the spin-orbit splitting), the excited state manifold splits into two triplets, each 
with a particular well defined spatial wavefunction. Orbital and spin degrees of freedom 
separate in this regime and spin-preserving transitions are excited by linearly polarized 
light. The strain Hamiltonian is given by ( [Rogers et al.^ 2009 Maze et a/., 2011) 


Hstrain — ^2 (|6a:)(^y| |^y)(^a:|) 


(1.9) 


where 61 and 62 are different parameters describing the crystal strain. Figure |1.3| 
shows an example of how 61 strain splits and mixes the excited states, and as a results 
changes the polarization properties of the optical transitions. A similar effect occurs 
for 62 strain, where IA 2 ) is mixed with |^i). 

The full Hamiltonian for the excited state manifold of the NV center is then 


H Hgg Hgo Hgixaini 


(1.10) 
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Fig. 1.3 Effect of strain on the properties of the NV center, a. Energies of the excited states 
as a function of strain, expressed in units of the linear strain induced splitting between the \Ex) 
and \Ey) states, b. Strain induced mixing of the \A2) and IE2) states, showing the fraction 
of A2 and E2 character of the highest energy excited state, c. Polarization character of the 
optical transition between IA 2 ) state and the ground state |+1) as a function of strain. As 
strain increases, the polarization changes from circularly to linearly polarized. The transition 
between IA 2 ) and | —1) shows the same behavior, except with cr_ and cr+ switched. Therefore, 
at high strain, decay from the IA 2 ) results in a separable state of the photon polarization and 
spin rather than an entangled state. 


which can be written in the set of basis states {|Ai), IA 2 ), , |E^i), |E^ 2 )} as 
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It is worth mentioning that while the excited state configuration is highly affected 
by strain, the ground state configuration is unaffected to first order due to its anti¬ 
symmetric combination of Cx and Cy orbitals. It is also protected by the large optical 
gap between the ground and excited state to second order perturbation in strain. 


1.3 Experimental techniques 

Having developed a theoretical model for the optical transitions of the NV center, we 
now move on to experimental explorations of these properties. While there are many 
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available descriptions of setups and procedures for room temperature experiments with 


NV centers (Gruber et a/., 1997 Childress et a/., 2006), we will focus here on some 


techniques relevant for low-temperature studies of the optical properties. 


1.3.1 Diamond materials 

We begin this section with a description of different types of single crystal diamond 
samples and their properties with regard to NV center properties. The first type of 
commonly available diamond material is classified as type Ib. In these diamonds, sub¬ 
stitutional nitrogen atoms (also called PI centers) are the dominant defects, which 
cause the samples to have a yellow color ( [Zaitsev, 2010 ). Type Ib diamonds can oc¬ 
cur naturally and are typically grown by the high-pressure-high-temperature (HPHT) 
method by companies such as Element Six and Sumitomo. Element Six specifies their 
type Ib diamond as having substitutional nitrogen concentrations of <200 ppm. In our 
experience, the concentration of NV centers in type Ib samples vary widely through 
different samples, but NV centers in bulk samples are usually unresolvable under con- 
focal microscopy. There have been no reports of coherent optical transitions from type 
Ib diamonds, except in one nanocrystal sample where a narrow, but spectrally unsta¬ 


ble line was observed under photoluminescence excitation (PLE) spectroscopy (Shen 
et a/., 2008). However, as far as we know, this result has never been reproduced. In 


our experiments, type Ib samples are typically used for developing and fabricating 
nanophotonic structures and studies of NV centers in these structures at room tem¬ 
perature. 

The rest of the diamond samples that are suitable for the experiments described 
here can be classified as type Ila, which is typically defined as having nitrogen impurity 
concentrations of less than Ippm. Type Ila diamonds are rare in nature, but one such 
sample was found containing a resolvable concentration of single NV centers that show 


optical linewidths that are close to the lifetime limited linewidth of 13 MHz (Tamarat 


et al.^ 2006). Many initial experiments on the optical properties of NV centers were 


performed on several different pieces of this sample, including the ones described in 


Sections 1.6 and 1.7 of this chapter, along with, for example, (Tamarat et a/., 2006) 


and (Batalov et a/., 2008). The unique properties of this sample, and its somewhat 
exotic origins in the Ural mountains, has earned it the nickname “The Magic Russian 


Diamond”. The piece used in our experiments is shown in Eignre pMp . Clearly, such 
specialized samples cannot be used for wide-scale experiments and development of 
diamond-based devices. An alternative is to work with synthetic type Ha diamonds 
grown using microwave assisted chemical vapor deposition (CVD) at companies such 
as Element Six Inc. Standard type Ha grade samples typically also show a unresolv¬ 
able concentration of NV centers. However, “electronic grade” diamonds developed by 
Element Six, which are specified to have substitutional nitrogen concentrations of <5 
ppb, have NV center concentrations that range from barely resolvable single NVs to 
less than I NV per tens of pm^. Native NV centers found in these samples often have 


linewidths that are broadened by spectral diffusion to several hundred MHz (Bernien 


et a/., 2012 Acosta et a/., 2012 Stacey et a/., 2012). However, they have been used in a 
variety of experiments involving optical transitions of NV centers (Bernien et al.^ 2013 
Bernien et a/., 2012 Sipahigil et a/., 2012), and are also excellent starting materials 
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Fig. 1.4 Experimental setup and characterization of NV centers, a. The type Ila natural 
diamond sample used for the experiments in Sections |1.6| and |1.7| b. Confocal microscope 
setup. Components in black are needed for a basic room temperature apparatus for investi¬ 
gating the ground state spin physics of the NV center. Red components, such as a cryostat 
and resonant excitation laser, are additionally needed for low temperature studies of coherent 
optical properties. A supercontinuum laser is added for characterization of photonic devices 
such as cavities and waveguides, and an additional, independent collection channel is used 
for collecting light from a different spatial location than the excitation, as shown in green. In 
addition, the spectrometer can also be used for photoluminescence studies for the NV center. 

c. Confocal image of single NV centers obtained using 532 nm excitation and PSB collection. 

d. Low temperature photoluminescence spectrum of a NV center. Sharp peak is ZPL at 637 


for introducing optically coherent NV centers through ion implantation, as described 
in Section fLSl 

1.3.2 Experimental setup for optical characterization of NV centers 

Figure |1.4|3 shows a schematic of a multi-purpose setup used for optical investiga¬ 
tion of NV centers and diamond-based photonic devices. Depending on the type of 
experiment, only parts of this comprehensive setup are needed. A typical basic con¬ 
focal microscope for room-temperature studies of NV centers is shown in black. A 
high numerical aperture (NA) objective is used to focus laser light onto the sample. 
In order to scan over the sample and locate individual centers in a confocal image 
such as shown in Figure |1.4| :, the diamond is either mounted on a high resolution 
positioning stage, or the optical path is scanned using a galvo mirror imaged onto 
the back of the objective with a pair of lenses comprising a 4f imaging system. The 
vibronic sidebands of the NV center’s optically excited states can be addressed using 
a laser of higher energy than the 637 nm zero-phonon line (ZPL). Typically, this is 
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done using a 532 nm frequency doubled NdiYAG laser, which is the most common 
method for controlling the charge state of the NV center. CW 532 nm excitation can 
ionize both the negatively charged NV center and surrounding charge donors such 
as substitutional nitrogens, resulting in a equilibrium charg e state of the NV center 
that is ^ 70% NV“ and ^ 30% NV^ (Aslam et al, 2013). Off resonant excitation 
is also used in room temperature experiments to polarize the spin of the NV center 
in to the \ms = 0) ground state. This process happens through non-spin preserving 
transitions from the optically excited state into spin-singlet states of the NV center, 
which then preferentially decays to the \ms = 0) state (Manson et a/., 2006). Finally, 
microwave and RF fields for driving spin transitions of the NV centers can be applied 
with a variety of methods, including copper wires, coplanar waveguides, and resonator 
structures. 

For low temperature experiments, the diamond sample is placed inside a helium 
flow cryostat (Janis ST-500) and cooled to ^4K. Alternative low temperature systems 
such as closed-cycle pulse tube coolers and bath cryostats can also used. Diamond 
samples are typically mounted onto the cold finger, or some intermediate substrate 
such as a silicon wafer, using indium or silver paste. For resonant excitation of the 
NV centers, a external cavity diode laser around 637 nm is introduced into the optical 
path. The New Focus TLB 6304 laser is a popular choice due to its large mode-hop- 
free scanning range of 60-80 GHz, wavelength tunability over several nanometers, and 
ease of operation. Both the 532 nm and 637 nm lasers need to be pulsed for most 
experiments, which can be accomplished using acoustic-optical modulators (AOM) or 
electro-optical modulators (EOM). 

Fluorescence from the NV center is collected through the same optical path, cou¬ 
pled into single-mode fibers, and detected using avalanche photodiodes (APD) or a 
spectrometer. Filters are placed in the collection path to remove reflected excitation 
light and separate different parts of the emission spectrum. As shown in Figure [L^, 
the NV center spectrum consists of a broad separate phonon-sideband (PSB) and a 
sharp ZPL at low temperatures. The percentage of emission into the ZPL is typically 
3 — 5%, which is a major limitation on the efficiency of applications based on coherent 
photons from the NV center. In Section |1.9[ we will address this issue by enhancing 
the emission into the ZPL using cavity QED with nanophotonic structures. In order to 
characterize these structures, we use an additional supercontinuum laser as a broad¬ 
band light source to perform transmission measurements of cavities and waveguides. 
The light is coupled into one port of the optical device, and the transmitted light must 
be collected from an output port at a spatially separate location. Eor this purpose, a 
second collection channel can be added with its own set of galvo mirrors. In addition, 
the signal from the two channels can be combined on a time-correlated single pho¬ 
ton counting (TGSPC) device for measurements of the second-order autocorrelation 
function g‘^{r) of the emission from individual NV centers to show that they are single¬ 
photon emitters (Kurtsiefer et a/., 2000). This can also be accomplished by splitting 
the light collected in a single channel using a fiber beamsplitter. 
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1.4 PLE spectroscopy of NV centers 

The optical transitions of the NV center can be identified using PLE spectroscopy by 
scanning the resonant excitation laser across the ZPL and collecting photons in the 
PSB. Since the |0) ^ \Ex) and |0) ^ \Ey) are best the cycling transitions, they are 
usually the ones observed in PLE spectroscopy in the absence of microwave driving or 
laser modulation. However, repeated resonant excitation eventually ionizes the NV“ 
and optically pumps the NV center into a different ground state spin sublevel due to 
the small non-spin conserving cross transitions out of the \ms =0) excited states into 
the |±1) ground states (Tamarat et a/., 2008). Therefore, some form of charge and spin 
repumping must eventually be performed. As shown in Eigure[L^, the first method 
we use involves scanning the 637 nm laser across the NV transitions and turning on 
the laser while simultaneous collecting PSB counts for ^ 2 ms at each frequency. 
At the end of the scan, the 637 nm laser is turned off and a 532 nm repumping 
pulse is applied. Any duration of repumping longer than a few /rs is sufficient to 
reinitialize the NV center. Depending on the charge environment of the NV center and 
its strain-dependent level structure, however, it may be difficult to obtain good signal 
from this simple method. Eor example, in the presence of a large number of other 
charge traps in the environment, the NV center may be ionized before a sufficient 
number of photons are scattered on resonance. Another potential issue is that the 
|0) ^ \Ex) and |0) ^ \Ey) transitions have some branching ratio into the rris = ±1 
states, causing a decrease in fluorescence on the timescale of several microseconds. 
However, fluorescence can potentially still be detected due to off-resonant excitation 
of transitions from the |±1) spin states. Eor example, at high strain, the |±1) ^ \Ai) 
transition becomes almost degenerate with the |±1) ^ \Ex) transition, allowing a 
steady state fluorescence to be established. However, this may not be the case for all 
NV centers. More robust methods of obtaining the PLE spectrum without repumping 
during the scan involve either CW microwave mixing of the |0) and |±1), or modulation 
of the resonant excitation laser by the ground state zero-field splitting, which also 
allows the \ms = ±1) transitions to be observed (Tamarat et a/., '2QQ8l|En et al, 2QQ9| 
Santori et a/., 200^ 


Eignre |1.5| 3 shows the results of PLE scans taken using this first method. A clear 
resonance can be seen in each scan, whose frequency then jumps during each succes¬ 
sive scan. This is an indication of spectral diffusion, where changes in the local electric 
field environment of the NV center causes the energies of its excited states to fluctuate 
over time. As mentioned earlier, from the electronic orbitals \ex) and \ey) involved 
in the excited states, one can see that the \Ex) and \Ey) states have a permanent 
electric dipole moment. This dipole moment has been measured to be ^0.8 Debye, 
which gives rise to energy shifts of ^4 GHz/MV/m (Maze et a/., 2011). It is now well 
established that 532 nm repumping causes complex charge dynamics in the diamond 


environment by photoionizing impurities that act as charge donors (Bassett et a/., 
2011). Therefore, while very narrow single-scan linewidths can be obtained using this 


method (Shen et a/., 2008 Tamarat et a/., 2006 Eu et a/., 2009), the long-term, ex¬ 


trinsically broadened linewidth obtained by averaging many such scans can be much 


broader (Bernien et a/., 2012 Acosta et a/., 2012 Stacey et al.^ 2012). This issue is 


addressed in more detail in Section p78] of this chapter. 
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Fig. 1.5 Methods for PLE spectroscopy of NV centers, a. Pulse sequence for repumping 
during each PLE scan (top), along with an example showing the result of successive scans 
(bottom) and the averaged spectrum (middle), b. Pulse sequence for repumping at the end 
of each PLE scan (top), along with corresponding results for the same transition as in a., 
except zoomed in over a smaller frequency range. 


A second method of obtaining the PLE spectrum involves a pulse sequence consist¬ 
ing of a 1 yus 532 nm repumping pulse, followed by a 10 /is resonant excitation pulse, 
during which PSB counts are collected. As the 637 nm laser is scanned across the NV 
transitions, this sequence is repeated many times for each frequency, as illustrated in 
Figure |1.5| This method has the advantage that it is robust against ionization and 
spin polarization, and should in principle give good signal to noise regardless of the 
charge environment and strain-dependent level structure of the NV center. It also 
shows the overall lineshape of the transition, including effects of spectral diffusion, in 
a single scan. Again, microwave fields are used to obtain a spectrum of all the possible 
transitions, including those that involve the = ±1) states. 


1.5 Optical properties of NV centers 


Figure [L^ shows a PLE spectrum taken with CW microwave excitation that mixes 
the ground states during the 637 nm laser pulse. All resonances expected from the 
direct and cross transitions are observed, as shown in |1.6| 3, along with some possible 


microwave-optical two photon transitions. 

The strongest transitions in 1.6r are the |0) ^ \Ex) and |0) ^ \Ey) cycling tran¬ 
sitions. Figure shows an example of optical Rabi oscillations on the |0) ^ \Ey) 
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Fig. 1.6 Optical transitions of the NV center, a. All possible transitions between the ground 
and excited states, with direct transitions indicated with solid lines, and spin non-conserving 
cross transitions indicated with dashed lines, b. PLE spectrum taken with CW microwave 
excitation, c. Frequencies of all possible transitions shown in a. as a function of strain. The 
frequencies of the peaks in b. are matched to a particular strain value. The extra unidentihed 
peak may be due to a two photon transition from the |0) state to the \Ex) state through the 
absorption of an optical photon and emission of a microwave photon. 

transition, measured by turning on a 40 ns resonant laser pulse a rise time of 
ns using an EOM (Guided Color Technologies). We then accumulate a histogram of 
photon arrival times on a TCSPC device with a timing resolution of 195 ps. The 
photon count rate is proportional to the population in the \Ey) state, which, in the 
absence of decay or decoherence, is given by Psyi^) =cos^(|^4|t), where ft = jl • E/fi 
is the Rabi frequency. In this case, however, the oscillations are damped by the finite 
lifetime of the excited states, which can be measured by fitting to the exponentially 
decaying tail after the pulse turns off. In addition to population decay, additional de¬ 
coherence, such as frequency fluctuations of the transition frequency relative to the 
laser frequency caused by, for example, spectral diffusion, will also contribute to the 
damping of Rabi oscillations. We note here, however, that this additional dephasing 
alone does not always give a good estimate of the spectral diffusion of the NV center, 
unlike what is sometimes claimed in literature. For example, if the transition frequency 
jumps well outside of the power-broadened linewidth of the transition, then the NV 
center will effectively not be excited, and there will simply be no contribution to the 
Rabi oscillation data. 
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Fig. 1.7 a. Optical Rabi oscillations detected using PSB fluorescence and resonant excita¬ 
tion, along with corresponding pulse sequence. The 637.2 nm laser is tuned to the |0) ^ \Ey) 
transition and is on from 50 ns to 90 ns. b. Rabi oscillations between |0) and | + 1) (blue) 
detected using resonant spin readout and corresponding pulse sequence. Red lines are cali¬ 
bration levels for the |0) and | + 1) states. To initialize the spin states for this calibration, we 
use a measurement-based technique where the |0) is prepared conditioned on the detection 
of a PSB photon after excitation to the \Ey) state. To prepare the | + 1) state, we then apply 
a microwave tt pulse before readout. 


By measuring Rabi oscillations, we can determine the length of laser pulse needed 
to, for example, perform an optical tt pulse to transfer population to the 1^42) state 
for entanglement generation, as described in Section [T6| Note that, while optical Rabi 
oscillations can also be observed between the \ms = il) ground and excited states, the 
dynamics are more complicated due to the fact that these are A-type transitions. In the 
absence of a magnetic field, instead of settling to a steady state value corresponding 
to 50% population in the excited state, the photon count rate continually decreases 
after a few optical cycles as population is pumped into a ’’dark” ground state that is 
not excited by the by the laser. However, the observed Rabi oscillations are usually 
sufficient for determining the optical tt pulse length. 

The cycling transitions of the NV center can also be used to efficiently read out 
the spin state of the NV center. As shown in Figure when a resonant laser pulse is 
turned on for 10 /is, the number of photons scattered is proportional to the population 
in the \ms = 0) states, and microwave-driven Rabi oscillations between the |0) and 
1+1) ground state sublevels can be observed. A detailed description of the resonant 
spin-readout procedure used in Section 1.6 can be found in (Togan et a/., “Mot - The 
efficiency and fidelity of this spin readout method is limited by the photon collection 
efficiency, the temperature, and the branching ratio out of the \ms =0) manifold. As 
can be predicted by the Hamiltonian given in Equation |1.2.3[ this branching ratio 
depends on the local strain. In our first experiments, the spin readout contrast was 
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limited to cm = 0.11 ± 0.0022 counts/shot when the spin was prepared in the |0) 
state and a background-limited cb = 0.0057 ± 0.0010 counts/shot when the spin 
was prepared in the |±1) states. Therefore, to determine the spin of the NV center, 
the state had to be prepared and read-out many times to obtain an average number 
of counts per shot C, which then allows us to calculate the population in the |0) 
state as P = (C — Cb) / (cm — cb). By increasing the collection efficiency with a solid 
immersion lens and working with low strain NV centers, single shot resonant readout 


of the spin state has subsequently been demonstrated (Robledo et a/., 2011). 


The rest of the chapter gives several examples of applications based on the optical 
properties of NV centers that we have just described. 


1.6 Quantum entanglement between an optical photon and a 
solid-state spin qubit 

1.6.1 Introduction 


A quantum network (Kimble, 2008) consists of several nodes, each containing a long- 


lived quantum memory and a small quantum processor, that are connected via entan¬ 
glement. Its potential applications include long-distance quantum communication and 
distributed quantum computation (Duan and Monroe, 2008). Several recent experi¬ 


ments demonstrated on-chip entanglement of solid-state qubits separated by nanome¬ 
ter (Neumann et a/., 2008) to millimeter (Ansmann et a/., 2009 DiCarlo et a/., 2009) 


length scales. However, realization of long-distance entanglement based on solid-state 


systems coupled to single optical photons (de Riedmatten et a/., 2008) is an outstand¬ 


ing challenge. The NV center is a promising candidate for implementing a quantum 
node. The ground state of the negatively charged NV center is an electronic spin triplet 
with a 2.88 GHz zero-field splitting between the \ms = 0) and = ±1) states (from 
here on denoted |0) and |±1)). With long coherence times (Balasubramanian et a/., 
2009), fast microwave manipulation, and optical preparation and detection ( [Fuchs 


et a/., 2009), the NV electronic spin presents a promising qubit candidate. Moreover, 


it can be coupled to nearby nuclear spins that provide exceptional quantum memories 
and allows for the robust implementation of few-qubit quantum registers (Neumann 


et a/., 2008 Dutt et a/., 2007). In this section we demonstrate the preparation of 


quantum entangled states between a single photon and the electronic spin of a NV 
center: 

|v|/) = T(|^_)|+l) + |^^)|_l)), (1.12) 

where |cr+) and |cr_) are orthogonal circularly polarized single photon states. 


The scheme we use here is inspired by previous work done in trapped ions (Bli¬ 


nov et a/., 2004). The key idea of our experiment is illustrated in Figure 1.8 The 


NV center is prepared in a specific excited state (IA 2 ) in Fig |1.8[ ) that decays with 
equal probability into two different long lived spin states (|±1)) by the emission of 
orthogonally polarized optical photons at 637 nm. The entangled state given by Eq. 
|1.12| is created because photon polarization is uniquely correlated with the final spin 
state. This entanglement is verified by spin state measurement using a cycling optical 
transition following the detection of a 637 nm photon of chosen polarization. 
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Fig. 1.8 Scheme for spin-photon entanglement. Following selective excitation to the \A 2 ) 
state, the A system decays to two different spin states through the emission of orthogonally 
polarized photons, resulting in spin-photon entanglement. 


1.6.2 Experimental demonstration of spin-photon entanglement 


We now turn to the experimental demonstration of spin-photon entanglement. The 


experimental setup is outlined in (Togan et a/., 2010). To create the entangled state, we 
use coherent emission within the narrow-band zero phonon line (ZPL), which includes 
only 4% of the NV center’s total emission. The remaining optical radiation occurs in the 
frequency shifted phonon side band (PSB), which is accompanied by phonon emission 
that deteriorates the spin-photon entanglement (Kaiser et a/., 2009). Isolating the weak 


ZPL emission presents a significant experimental challenge due to strong reflections of 
the resonant excitation pulse reaching the detector. By exciting the NV center with a 
circularly polarized 2 ns tt pulse that is shorter than the emission timescale, we can 
use detection timing to separate reflection from fluorescence photons. A combination 
of confocal rejection, modulators, and finite transmittivity of our optics suppresses the 
reflections sufficiently to clearly detect the NV center’s ZPL emission in a 20 ns region. 

For photon state determination, ZPL photons in either the \a±) or \H) = ^(|^+) + 
|cr_)), \V) = ^(|cr+) — |cr_)) basis are selected by a polarization analysis stage and 
detected after an optical path of ^ 2 m. Spin readout then occurs after a 0.5 /is spin 
memory interval following photon detection by transferring population from either 
the |±1) states or from their appropriately chosen superposition into the |0) state 
using microwave pulses fl±i. The pulses selectively address the |0) |±1) transitions 

with resonant frequencies uj± that differ by Suj = — uj- =122 MHz due to an 

applied magnetic field. For superpositions of |±1) states an echo sequence is applied 
before the state transfer to extend the spin coherence time. The transfer is followed by 
resonant excitation of the |0) \Ey) transition and collection of the PSB fluorescence. 
We carefully calibrate the transferred population measured in the |0) state using the 
procedure detailed in ( [Togan et a/., 2QTQ ). Figure 1.9r shows the populations in the 
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|±1) states, measured conditionally on the detection of a single circularly polarized 
ZPL photon. Excellent correlations between the photon polarization and NV spin 
states are observed. 





Fig. 1.9 Measurement of spin-photon correlations in two bases, a. Conditional probability 
of measuring |dil) after the detection of a cr-h or a— photon, b. Conditional probability of 
measuring |zb) after the detection of a if or E photon, extracted from a fit to data shown 
in (c) and (d). c&;d. Measured conditional probability of finding the electronic spin in the 
state \M) after detection of a E (c) or H (d) photon at time td- Blue shaded region is the 
68% confidence interval for the fit (solid line) to the time-binned data. Errors on data points 
are one STD. Combined with the data shown in Figure 1.9'a), oscillations with amplitude 


outside of the yellow regions result in fidelities greater than 0.5. The visibility of the measured 
oscillations are 0.59 ± 0.18 (c) and 0.60 ± 0.11 (d). Reprinted from (Togan et al, 2010). 


To complete the verification of entanglement, we now show that correlations persist 
when ZPL photons are detected in a rotated polarization basis. Upon detection of a 
linearly polarized \H) or |E) photon at time td^ the entangled state in Eq 1.12| is 


projected to |±) = ^ (|+1) i l“l ))5 respectively. These states subsequently evolve in 
time (t) according to 


|±)^ = T 1+1) ±e-“-(*-*<*) 1-1)) . 


(1.13) 


In order to read out the relative phase of superposition states between |+1) and | —1), 
we use two resonant microwave fields with frequencies and uu- to coherently transfer 
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the state |M) = ^ e-h^-^-(0+-0-)) |_i^^ to |0), where the initial 

relative phase 0+ — is set to the same value for each round of the experiment. 
Thus, the conditional probability of measuring the state |M) is 


PM\H,V (td) — 


1 ± COS a (td) 


(1.14) 


where a (td) = (ci;+ — oo-) td + (0+ — 0_). Eq. 1.14 indicates that the two conditional 
probabilities should oscillate with a tt phase difference as a function of the photon 
detection time td. This can be understood as follows. In the presence of Zeeman split¬ 
ting {6uj 7 ^ 0) the NV center’s spin state is entangled with both the polarization and 
frequency of the emitted photon. The photon’s frequency provides which-path informa¬ 
tion about its decay. In the spirit of quantum eraser techniques, the detection of \H) or 
\V) at td with high time resolution 300 ps <C l/5uj) erases the frequency information 
( [Scully and Drnhl, 1982 Wilk et al.^ 2007). When the initial relative phase between 
the microwave fields ft±i is kept constant, the acquired phase difference (cc;+ — uj-)td 
gives rise to oscillations in the conditional probability and produces an effect equiva¬ 
lent to varying the relative phase in the measured superposition, allowing us to verify 
the coherence of the spin-photon entangled state. 

The detection times of ZPL photons are recorded during the experiment without 
any time gating, which allows us to study spin-photon correlations without reducing 
the count rate. The resulting data are analyzed in two different ways. First, we time- 
bin the data and use it to evaluate the conditional probabilities of measuring spin state 
\M) as a function of \H) or \V) photon detection time (Figure jL^ c,d). Off-diagonal 
elements of the spin-photon density matrix are evaluated from a simultaneous fit to 
the binned data and the time bins are chosen to minimize fit uncertainty as described 


in (Togan et al.^ 2010). The resulting conditional probabilities are used to evaluate a 


lower bound on the entanglement fidelity of F > 0.69 ±0.068, above the classical limit 
of 0.5, indicating the preparation of an entangled state. 

We further reinforce our analysis using the method of maximum likelihood esti¬ 


mate. As described in (Togan et a/., 2010), this method is applied to raw, un-binned 


ZPL photon detection and spin measurement data and yields a probability distribu¬ 
tion of a lower bound on the fidelity. Consistent with the time-binned approach, we 
find that our data are described by a near Gaussian probability distribution associated 
with a fidelity of F > 0.70 ± 0.070. Significantly, the cumulative probability distribu¬ 
tion directly shows that the measured lower bound on the fidelity is above the classical 
limit with a probability of 99.7 %. 

Several experimental imperfections reduce the observed entanglement fidelity. First, 
the measured strain and magnetic field slightly mixes IA 2 ) state with the other ex¬ 
cited states. We estimate that IA 2 ) state imperfection and photon depolarization in 
the setup together reduce the fidelity by 12 %, the latter being the dominant effect. Im¬ 
perfections in readout and echo microwave pulses decrease the fidelity by 3 %. Other 
error sources include finite signal to noise in the ZPL channel (fidelity decrease 11 
%), as well as timing jitter (another 4 %). The resulting expected fidelity (73 %) is 
consistent with our experimental observations. Finally, the entanglement generation 
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succeeds with probability p ^ 10“^, which is limited by low collection and detection 
efficiency as well as the small probability of ZPL emission. 

The demonstration of spin-photon entanglement is the first step in implementing 
a quantum network architecture for long-distance quantum communications using NV 


centers (Kimble, 2008). Such an entangled state connects the solid state, spin-based 


local quantum register with optical photons, which can then be disseminated over free 
space or through optical fibers. Following the above experiment, tremendous progress 
has been made in achieving the next steps in building a diamond-based quantum net¬ 
work. Hong-Ou-Mandel interference was demonstrated with indistinguishable photons 


from remote NV centers hosted in separate diamond samples (Sipahigil et a/., 2012 
Bernien et a/., 2012). Combined with spin-photon entanglement, this recently allowed 


for the remote entanglement of two NV centers separated by a distance of three meters 


through entanglement swapping (Bernien et a/., 2013). 

1.7 Laser cooling and real-time measurement of nuclear spin 
environment of a solid-state qubit 

In this section, we describe the use of the NV center’s optical transitions and coherent 
population trapping (CPT) to sense and manipulate its surrounding spin bath, which 
consists of the spin of the ^^N nucleus associated with the NV center itself and the 
nuclei in the diamond lattice. Over the past two decades, CPT has been employed 


for laser cooling of neutral atoms and ions ( 

Aspect et a/., 1988), creation of ultra- 

cold molecules (Ni et al.^ 2008), optical magnetometry (Scully and Fleischhauer, 1992 

Budker and Romalis, 2007 

), and atomic clocks (Vanier, 2005), as well as for slowing 

and stopping light pulses ( 

Fleischhauer et aL, 2005 

). The electronic spin of the NV 


center is a promising system for extending these techniques to the solid state. The 


NV center has a long-lived spin triplet as its electronic ground state (Manson et aL, 


2006), whose iris = ±1,0 sublevels are denoted as |±1) and |0). In pure samples, the 


electron spin dynamics are governed by interactions with the spin-1 ^^N nucleus of 
the NV center and spin-1/2 ^^C nuclei present in 1.1 % natural abundance in the 
diamond lattice (Fig. |1.1Q^ ). Control over nuclear spins (Dutt et aL, 2007 Neumann 
et g/., 2010 ) is of interest for both fundamental studies and for applications such as 
nanoscale magnetic sensing ( [Maze et al.^ 2008 Balasubramanian et a/., 2008) and 
realization of quantum networks (Childress et a/., 2006 Togan et a/., 2010). Here we 
achieve such control via two complimentary methods: effective cooling of nuclear spins 


through nuclear state selective CPT (Issler et a/., 2010) and conditional preparation 


based on fast measurements of the nuclear environment and subsequent post-selection 


(Giedke et aL, 2006). 


While most prior work involved the use of microwave and RF fields for manipu¬ 
lating both the electronic and nuclear spin states, we utilize all-optical control of the 


electronic spin (Robledo et al.^ 2010 Santori et a/., 2006 Buckley et a/., 2010). Specif¬ 


ically, we make use of A-type level configurations involving the NV center’s \Ai) and 
IA 2 ) optically excited electronic states and the |±1) ground states (Fig. |l.lQ^ ) ( [Togan 
\et aL, 20TQ||Maze et aL, 2011 ). At low temperatures (<10 K) and in the limit of zero 
strain, \Ai) and | A 2 ) are entangled states of spin and orbital momentum coupled to the 
I-hi) (| —1)) state with cr_ (cr+) circularly polarized light. Correspondingly, excitation 
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with linearly polarized light drives the NV center into a so-called dark superposition 


state when the two-photon detuning is zero (Scully and Zubairy, 1997). In the present 


case the two-photon detuning is determined by the Zeeman splitting between the |±1) 
states due to the combined effect of the Overhauser field originating from the nuclear 
spin environment and any externally applied magnetic field (Issler et a/., 2010 [Stepa- 
nenko et a/., 2006). When the external field exactly compensates the Overhauser field, 


the electronic spin of the NV center is pumped into the dark state after a few optical 
cycles and remains in the dark state, resulting in vanishing fiuorescence. This is the 
essence of the dark resonances and CPT. 


1.7.1 Coherent population trapping with NV centers 

In our experiments, the \Ai) and | A 2 ) states are separated by approximately 3 GHz and 
are addressed individually with a single linearly polarized laser at near zero magnetic 
field. Since there is a finite branching ratio from the rus = ±1 manifold of the electron 
spin into the |0) state, we use a recycling laser that drives the transition between |0) 
and the \Ey) excited state, which decays w ith a small but non-vanishing probability 
10“^) back to the |±1) states. Figure 1.10 3 presents experimental observation 


( _ 

of the CPT spectrum as a function of an external magnetic field at three different 
powers of a laser tuned to the |±1) ^ IH 2 ) transition. While a broad resonance 
is observed at high power levels, as the power is reduced, we clearly resolve three 
features in the spectrum separated by 4.4 MHz, which is two times the hyperfine 
splitting between three nuclear spin states. This separation corresponds to the 
magnetic field required to bring the electronic rris = ±1 hyperfine states with equal 
nuclear spin projection {rui = ± 1 , 0 ) into two-photon resonance. 

The dependance of the CPT resonance width upon the laser power, shown in 
Figure [iT^, reveals an important role played by repumping on the near-cycling |0) 
\Ey) transition. In contrast to a conventional, closed three-level system, this recycling 
transition can be used to enhance the utility of our CPT system by both decreasing 
the width of the CPT resonance and increasing the signal to noise ratio. The \Ai) state 
decays into the rUs = 0 ground state through the singlet with a substantial probability 
of ^ 40%. However, the population is returned back into the = ±1 state from 
\Ey) only after ^ 100 optical excitation cycles. As a result, away from the two-photon 
resonance, the NV quickly decays to the |0) state after being excited, where it then 
scatters many photons through the |0) \Ey) cycling transition before returning to 
the A system. If the NV center is not in a dark state, this process effectively increases 
the number of photons we collect by 2 /r] = where 7 ce is the cross transition 

rate from \Ey) into |±1) and 7^1 is the rate from \Ai) to the singlet. The cycling 
effect also reduces the width of the CPT line since the |0) \Ey) transition is quickly 
saturated away from two-photon resonance, provided that the CPT laser excitation 
rate exceeds the leakage rate out of recycling transition. Significantly, both of these 
effects lead to improved sensitivity of dark resonances to small changes in two-photon 
detuning. 

To demonstrate this effect, the widths of dark resonances observed via excita¬ 


tion of |Ai) and IA 2 ) are compared in Figure 1.10Through an independent mea¬ 
surement of the branching ratios (see (Togan et a/., 2011)), we determined that 
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Fig. 1.10 Coherent population trapping in NV centers a. The A-type transitions between 
the ground states |±1) and excited states |Ai, 2 ) of a single NV center are addressed with a 
CPT laser, while a recycling laser drives the |0) to \Ey) transition. An external magnetic held 
is applied using a solenoid, b. Photon counts from NVa in a 300 ps window are plotted versus 
the applied held for 10 pW (blue), 3 pW (red), and 0.1 pW (yellow) of laser power addressing 
the IA 2 ) state. Blue and red datasets are shifted vertically by 5 and 2 counts for clarity, c. 
Width of individual CPT lines versus CPT laser power when the \Ai) (blue) or IA 2 ) 
(red) state is used. Errorbars in all hgures show ±1 s.d. Solid curves represent theoretical 
models. Reprinted from (Togan et al., 2011). 


1 ^ 1 ( 2 )) corresponds to an open (nearly closed) A system with ^ 3.1 x 10“^ 
{r]A 2 ^ 2.6). These experimental results are compared with a theoretical model de¬ 
scribed in ( [Togan et a/., 2 QlT ), which predicts that the resonance linewidth (5o is given 

by <5o = JRa/ 1 + ^ (^ + ^) PRaRetiI/{Re + 1 ) for small ??, where Ra(e) 

corresponds to the optical excitation rate by a laser tuned to the A{E) state and 7 
is the decay rate of excited states. The width at low powers is determined by the 
random magnetic held associated with surrounding ^^C nuclear states. When this line 
broadening mechanism is taken into account, the experimental results are in excellent 
agreement with these predictions, plotted as solid lines, showing that a high degree of 
optical control over the electronic spins of NV center can be achieved. 

Subsequent experiments (Togan et aL, 2011) showed that CPT can be used for 
not only probing the nuclear spin environment of the NV center, but also controlling 
and preparing the state of the spin bath. For the ^^N spins, it is possible to optically 
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cool the nuclear spin states using dark resonances. The physics is reminiscent of laser 


cooling of atomic motion via velocity-selective CPT (Aspect et a/., 1988 Issler et al.^ 


2010). A redistribution of the spin state population upon optical excitation takes 


place because the hyperfine coupling in the excited electronic state of the NV center 
is enhanced by a factor of ^ 20 compared to the ground state (Fuchs et a/., 2008). 
If the external field is set such that, for example, the rui = 0 hyperfine states are in 
two photon resonance, only the states with nuclear configuration mj = ±1 will be 
promoted to the excited states, where flip-flops with the electron spin will change the 
nuclear spin state to mj = 0. When the NV center spontaneously decays into the 
dark superposition of electronic spin states, optical excitation will cease, resulting in 
effective polarization (cooling) of nuclear spin into rui = 0 state. 



Fig. 1.11 Fast measurement and preparation of spin bath. a. Counts from 200 successive 
fast magnetic field ramps are shown on horizontal lines. The centers of Lorentzian fits to 
individual runs are indicated with a dot. b. Pulse sequence for preparation and subsequent 
measurement of the configuration. Bprep is set within the central line while Bjio is 
varied to cover all associated states. Counts during a conditioning window of length Tcond 
(^cond) at the end of preparation and the readout window (uro) are recorded for each run. 
Data presented is an average of many such experimental runs. b. uro versus Bjio is shown 
in the middle plot with double Lorentzian fit. The same dataset analyzed by keeping only 
events with Ucond = 0 is shown in the bottom plot. Unprepared ^^C distribution is shown in 
the top plot for comparison (shifted by 4.3 counts for clarity). Reprinted from (Togan et al, 
2011 ). 


We also extended this technique to control the many-body environment of the NV 
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center, consisting of nuclei distributed throughout the diamond lattice. The large 
number of nuclear spin configurations associated with an unpolarized environment re¬ 
sults in a random Overhasuer field (5ov) with unresolved hyperfine lines. It produces 
a finite CPT linewidt h in m easurements that average over all configurations of the 
spin bath (Figure 1.103, c). However, one can make use of fast measurements 
and the long correlation time associated with evolution of the nuclear bath to 

observe its instantaneous state and and manipulate its dynamics. For example, such 
a fast measurement is illustrated in Figure where the externally applied field 

is ramped across a single mi = 0 line while the CPT lasers are on, and counts 
collected in 80 ps time bins are plotted for successive individual runs, many of which 
distinctly show a narrow dark region. Lorentzian fits to selected experimental scans 
reveal “instantaneous” CPT resonances with linewidths that are over a factor three 
less than those of the averaged measurement. The motion of the dark line centers 
(green dots in Figure 1.11^ ) indicates that the instantaneous field evolves in time. 
Fast measurements can also be used to conditionally prepare the ^^C environment 
of the NV center in a desired state with post-selection. Here, the nuclear spins are 
first prepared though optical excitation at a fixed value of the external magnetic field 
^prep 5 as shown in Figure l.llj p. By conditionally selecting zero photon detection 
events during this preparation step, we can select the states of the ^^C environment 
with vanishing two photon detuning 6 = 2^/iB(^prep + ^ov) = 0, where is the Bohr 
magneton. The resulting modified distribution of Overhauser field is then probed in a 
readout step by rapidly the changing magnetic field value to which is scanned 

shows (unconditioned) read- 


1.11 


across the resonance. The middle curve in Figure 
out counts recorded following the preparation step, while the bottom curve shows 
the result of measurement-based preparation. The measured width of such a condi¬ 
tionally prepared distribution is significantly smaller than the width corresponding to 
individual resonances obtained without preparation. 


1.8 Coherent optical transitions in implanted nitrogen vacancy 
centers 

Since one-of-a-kind natural diamond samples such as the one investigated in the pre¬ 
vious sections cannot be used as reproducible starting materials for scalable systems, 
recent experiments have made use of NV centers incorporated into high quality syn¬ 


thetic diamond during the growth process (Bernien et a/., 2012 Acosta et a/., 2012 


Stacey et a/., 2012). In addition to imperfections in their optical coherence, NV centers 


introduced during the diamond growth process are limited in their usefulness due to 
their low concentration and random locations. For example, in order to incorporate 
NV centers into nanophotonic or nanomechanical devices, it is necessary to develop a 
method to create a controllable concentration of emitters in a well defined device layer 
while maintaining their excellent optical properties. 

This section describes how NV centers can be created in a surface layer that exhibit 
coherent optical transitions with nearly lifetime-limited linewidths even under the ef¬ 
fects of spectral diffusion. We accomplish this by first introducing NV centers at a well 
defined and controllable depth using ion implantation (Pezzagna et a/., 2010). We then 
demonstrate how to suppress the fundamental cause of spectral diffusion by creating 
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a diamond environment that is nearly free of defects that contribute to charge fluctu¬ 
ations. This is achieved through a combination of annealing and surface treatments. 
Finally, we further improve the optical properties of the NV centers by combining 
our approach with a newly developed technique for preventing photoionization of the 
remaining charge traps (Siyushev et a/., 2013). 
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Fig. 1.12 Procedure for creating spectrally stable NV centers through ion implantation, 
annealing, and surface treatment. The Ar/Ch (O 2 ) etch was done using an ICP power of 400 
W (700 W), a RF power of 250 W (100 W), a DC bias of 423 V (170 V), a chamber pressure 
of 8 mTorr (10 mTorr), and a gas flow rate of 25/40 seem (30 seem). The sample temperature 
was set to 17°C. Reprinted from (Chu et al. 2014). 


Fig. ra a) summarizes the experimental procedure to create shallow implanted 
NV centers with narrow optical lines. As our starting material we use single crystal 
diamond from Element Six grown using microwave assisted chemical vapor deposition 
(CVD). The single crystal samples were homoepitaxial grown on specially prepared 
(100) oriented synthetic diamond substrates, taking care to ensure that the surface 
quality of the substrates was as high as possible to reduce sources of dislocations. The 


CVD synthesis was performed using conditions as described in (Isberg et al.^ 2002 


Markham et al, 2011). The resultant diamond material is measured using electron 


paramagnetic resonance (EPR) to have a nitrogen concentration of less than 5 parts 
per billion. The surface layer of these samples is generally damaged and highly strained 
because of mechanical polishing. Therefore, we first remove the top several micrometers 
of the sample using a 30 min Ar/Cl 2 etch followed by a 20 min O 2 etch in an UN AXIS 
Shuttleline inductively coupled plasma-reactive ion etcher. We then implant ^^N+ ions 
at an energy of 85 keV at a variety of doses from 10^/cm^ to 10^^/cm^. According 
to simulations done using Stopping and Range of Ions in Matter (SRIM), this should 
result in a mean nitrogen stopping depth of ^100 nm with a straggle of ^20 nm 
(jZiegler, 2013). Subsequently, we anneal the implanted samples under high vacuum 
{P < 10 fcj Xorr). Our standard annealing recipe consists of a 4 hour 400°C step, 
a 2 hour 800°C step, and a 2 hour 1200°C step, with a 1 hour ramp up to each 
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temperature. The motivation behind these temperature steps will be explained later 
in the text. It is important that a good vacuum is maintained throughout the annealing 
process, since diamond etches and graphitizes under residual oxidizing gasses at such 
high temperatures (Seal, 1963). Following the anneal, we remove graphitic carbon 


and other surface contaminants by cleaning the diamonds for one hour in a 1:1:1 
refluxing mixture of sulfuric, nitric, and perchloric acids. Finally, we perform a 465°C 
anneal in an O 2 atmosphere in a rapid thermal processor (Modular Process Technology, 
RTP-600xp) in three 48 minute steps following the recipe used in (Fu et a/., 2010). 
This oxygen anneal is believed to remove sp‘^ hybridized carbon and result in a more 
perfect oxygen termination of the surface than acid cleaning alone. We have observed 
that it enhances charge stability in the NV“ state and helps to reduce blinking and 

photobleaching. _ 

Let us focus on samples implanted with a ion dose of 10^^/cm^. Fig. 1.12[ b) 
shows a histogram of linewidths for 50 NV centers in four such samples. While there is 
a range of linewidths, a majority of NV centers exhibit linewidths less than 500 MHz, 
which is comparable to naturally occurring NV centers deep inside electronic grade 


CVD diamond (Bernien et a/., 2012). This is an indication that we have eliminated 


most impurities with fluctuating charge states in the vicinity of the NV centers, which 
may include both impurities on the nearby surface and defects inside the diamond 
introduced by the implantation process. In Fig. 1.12[ c), we present an example spec¬ 
trum of the narrowest linewidth we have observed, along with the results of twenty-one 
successive PLE scans to demonstrate the long-term stability of the transition line. Our 
measured linewidth of 37 ± 4.8 MHz under conventional 532 nm repumping is, to the 
best of our knowledge, a record for the extrinsically broadened linewidth of NV centers. 


1.9 Diamond-based nanophotonic devices for quantum optics 

1.9.1 Cavity QED with NV centers: the idea 

In the previous sections, we have demonstrated the potential of using the NV center’s 
optical properties for a variety of different applications. However, the NV center is 
not a perfect optical emitter for several reasons. First, the Debye-Waller factor 
defined as the ratio of emission into the ZPL vs. the total emission, is only 3-5 %. 
This is dictated by the relatively large shift in nuclear coordinates between the NV 
center’s ground and excited states, which, by the Franck-Condon principle, results 
in decay from the excited state into phonon-coupled higher vibrational levels of the 


ground state (Gali et a/., 2011). In addition to such intrinsic properties of the NV 


center, there are also external effects such as spectral diffusion that result in non¬ 
ideal optical properties. While we can eliminate these effects to a large degree, as 
demonstrated in Chapter |1.8[ some residual imperfections in spectral stability usually 
remain. Additional techniques to stabilize the transition frequency, such as feedback 
techniques with applied electric fields, reduce the repetition rate of the experiments. 
Finally, it is generally difficult to extract photons from bulk diamond material with 
high efficiency. Due to the high index of refraction, total internal reflection confines 
most of the light inside the diamond at an air interface. In general, as shown in 
(Plakhotnik et a/., 1995), the collection efficiency of a lens or objective with numerical 
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Fig. 1.13 Collection efficiency of light from a dipole emitter, a. Collection efficiency vs. 
dipole orientation for NA=0.95 and n 2 =l. b. Collection efficiency vs. dipole orientation for 
NA=0.95 and n2=2.4. c. Collection efficiency vs. NA for 0em = 7r/2 and n 2 = 2.4. d. 
Collection efficiency vs. n 2 for 0em = 7r/2 and NA=0.95. 

aperture NA in a medium with index ni for light from a dipole source in index n 2 is 
given by 

C = i (4 - ScOsOmax “ COS^Omax + ^COS^Omax “ COsOmax)(t>em) , (I-IS) 

where 0em is the angle of the dipole axis with the optical axis, and Omax = sin“^(NA* 
^ 1 /^ 2 ) is the effective NA of the collection optics. For reference, the collection efficiency 
is plotted against several parameters in Figure [T.13| As can be seen, even for optimal 
dipole orientation and a 0.95 NA objective in air, only around 6% of the emitted 
photons are collected from a NV center in bulk diamond. 

The combination of these factors result in a very low rate of indistinguishable, co¬ 
herent photons from single NV centers. One promising path toward alleviating all of 
these imperfections is incorporating NV centers with photonic structures. For exam¬ 
ple, optical cavities can be used to enhance the emission into the ZPL and overcome 
spectral diffusion, while coupling to waveguide structures can increase the efficiency 
of photon collection and detection. 

The idea of using cavity QED to improve the emission properties of a NV center 
is as follows: When a NV center is placed inside the an optical cavity, its emission 
into the cavity mode is enhanced. By tuning the cavity mode to the NV center’s 
ZPL, the ZPL emission is therefore enhanced relative to the PSB emission. When 
this enhancement becomes larger than the Debye-Waller factor, the ZPL becomes 
the dominant emission channel of the NV center. This enhancement of the emission 
also means that the lifetime of the NV center is decreased and the lifetime-limited 
linewidth of the ZPL is increased. When the lifetime-limited linewidth becomes much 
larger than the spectral diffusion, the NV center can be viewed as a perfect emitter of 
indistinguishable photons. In other words, the energy uncertainty of the excited state 
caused by fast decay into the cavity mode is large enough that the energy uncertainty 
introduced by spectral jumps is negligible. 

We can calculate lifetime of a NV center coupled to the cavity (r) compared to 
one that is decoupled from the cavity (tq). They are related through a quantity called 
the Purcell factor, given by 



(1.16) 
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where Xo/n is the wavelength of the resonant mode inside the cavity and Q is the 
quality factor of the cavity. Ejyy and Ej^ax are the electric fields at the NV and 
maximum field inside the cavity, respectively. The mode volume of the cavity is defined 
as 

(\Eix)\^dV 

V=^-E-^ - (1.17 

\E 

Finally, pzPL is the effective dipole moment of the ZPL transition, which is the only 
part coupled to the cavity, while fiNv is the total dipole moment of the NV center. 
The Debye-Waller factor is then given by ^ = l/izPLp/lMAry P- The Purcell factor gives 
the enhancement of emission into the cavity mode over the emission when the NV is 
decoupled from the cavity. Here, we define “decoupled” as an NV that is still spatially 
located in the cavity, but detuned from the cavity mode in frequency, as is the case 
in our experiments. This means that n in the expression above is the effective index 
including all possible modes of the structure. Although this quantity can be difficult 
to estimate, it is somewhere between 1 and 2.4, the index of bulk diamond. 

While the Purcell factor expresses the relevant enhancement of emission in terms 
of experimentally relevant quantities, it is equivalent to the cooper at ivity, which gives 
more physical intuition as simply the ratio between the coupling strength of the emitter 
to the cavity and the decay rates of the system. The Purcell factor or cooperativity 
can also be written as 

P=— (1.18) 


where k = uj/Q is the decay rate from the cavity, with uj being the cavity frequency. 
7 is the decay rate of the NV center into modes other than the resonant cavity mode, 
and is given by 


^ Sirehc^ 3eh V^o/ 


(1.19) 


The single-photon Rabi frequency g is given by ^ = E^y • PzplI^- We emphasize here 
that since the cavity mode only couples to the ZPL, we use the effective dipole moment 
of that transition only. Inserting these expressions into Equation |1.18| and using the 
fact that e J \E{x)\‘^dV = Tiuo^ we recover the expression for the Purcell factor given in 
Equation |1.16| 

We note here that there are several alternative definitions of the Purcell factor 
that might also relevant here. One can also consider the lifetime change relative to 
emission into bulk diamond or free space, in which case the index n in Equation 


ZPL emission into the cavity relative to the ZPL emission when the NV is decoupled 
from the cavity. This is given by P' = P/^, and a useful quantity when comparing the 
Purcell enhancement to other cavity QED systems where the emitter has a single decay 
channel. In other words, it gives a better comparison for the quality of the emitter- 
cavity coupling without factoring in the internal imperfections of the emitter itself, 
and is the quantity usually quoted in the literature. It is also the relevant quantity for 
estimating the enhancement of the ZPL photon rate from the cavity compared to, say, 
bulk diamond. 


1.16 should be 2.4 or I, respectively. Another quantity is the enhancement of the 
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As explained above, in order to radiatively broaden the NV center through Purcell 
enhancement to overcome spectral diffusion, we require that P ^ or P' ^ 
where P is the extrinsically broadened linewidth. From Equation |1.16[ we see that to 
obtain a large Purcell factor, we need large Q and small mode volume. This can be 
achieved using photonic crystal cavities, where the mode volume can be designed to 
be ^ (A/n)^. Then, for a NV center with F ^ 100 MHz, a cavity with a Q of ^ 3000 
is required. We now discuss some directions for the physical implementation of such 
cavities in diamond. 


1.9.2 Diamond-based nanophotonic devices 

Over the past several decades, a large range of techniques have been developed for the 
fabrication of nanoscale devices in semiconductor materials ranging from silicon-based 


materials to III-V compounds such as gallium arsenide (England et a/., 2007; Blanco 


et a/., 2000). While many of these techniques can be transferred to the fabrication of 


diamond, one major difference between diamond and most semiconductor materials 
presents a significant challenge to the fabrication of nanoscale devices for confining 
light. Because of the large lattice mismatch of diamond with most other crystalline 
materials, it is not possible to grow high quality thin films of single crystal diamond 
on top of a different substrate. This means that we cannot take advantage of the 
index mismatch between the device and substrate to confine light through total inter¬ 
nal reflection. Here, we will briefly describe two methods for creating diamond-based 
nanophotonic devices by either creating a thin film that is placed on a substrate, or 
directly fabricating suspended structures on the surface of a bulk diamond sample. 
The first method for fabricating diamond-based photonic crystal cavities begins with 
a mechanically polished diamond membrane that is thinned down using a deep re¬ 
active ion etch (RIE) in inductively coupled plasma (ICP) of oxygen to the device 
thickness (Hausmann et a/., 2013 Earaon et a/., 2QTT] ). Devices are then defined on 
these thin films by e-beam lithography (EBL) using hydrogen silsesquioxane (HSQ) 
resist, which forms a mask that is transferred into the diamond film during another 
RIE etching step. The second technique for scalable fabrication of nanophotonic struc¬ 
tures involves undercut structures suspended above the surface of bulk diamond. This 


is achieved using a technique described in (Burek et al.^ 2012), where the diamond is 


placed inside a faraday cage during RIE so that the ions are redirected and impinge 
on the surface at an angle. As a result, instead of etching the diamond in the direction 
perpendicular to the surface, it undercuts the material beneath the mask from two 
directions, resulting in suspended structures with a triangular cross-section that are 


suspended from or supported by anchors. As shown in (Burek et al.^ 2012) and Eigure 


|1.14[ we can make a variety of devices including ring resonators, cantilevers, and ID 
photonic crystal cavities using this method. In addition, ID waveguide structures can 
be incorporated into hybrid cavities where a Bragg structure is defined by another 
material around an NV inside the waveguide, such as PMMA. 


1.9.3 Example: Purcell enhancement from a coupled NV-cavity system 

As an example, we present here a measurement demonstrating the Purcell enhance¬ 
ment of a NV center’s ZPL emission by a nanoscale hybrid diamond-PMMA cavity 
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Fig. 1.14 Suspended diamond-based nanophotonic structures, a. Procedure for angled RIE 
etching. A mask is defined using EBL on the surface of the diamond. A top-down etching step 
performed, followed by an angle-etching step where the diamond is placed in the center of a 
machined aluminum faraday cage shown on the bottom left. The mask is removed, leaving a 
suspended structure attached at both ends to anchors (not shown) b. SEM image of a large 
array of ~1000 angle-etched waveguides, c-d. SEM images of photonic crystal cavities formed 
in ID waveguides, e. Transmission spectrum of one such device. See (Burek et al, 2012) for 
details. 
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(de Leon et al.^ 2014). These measurements were made on a NV center with a 3 GHz 
linewidth inside a hybrid cavity with a Q of ^600. 

In order to tune the cavity mode frequency, we use a combination of two techniques. 
The first method uses a O 2 plasma etcher (Technics, Model 220) to remove a layer of 
material from the structure ( Hausmann et al.^ 2013| ). This tunes the cavity resonance 
to lower wavelengths because more of the cavity field resides in air than before and 
therefore has a higher energy. We have been able to achieve a wavelength shift of ^15 
nm without significant degradation in the Q factor using this method. Complementary 
to this blue tuning process, we can also red tune the cavity resonance by depositing 


neon gas on the structure when it has been cooled down inside the cryostat (Mosor 


et a/., 2005). We note here that gas deposition is a common method for tuning photonic 


crystal cavities, and a tuning range of ^23 nm was observed in (Hausmann et al.^ 2013). 


However, it is somewhat less effective for our hybrid structures since the high aspect 
ratio of the PMMA slabs prevent the diffusion of gas onto the diamond waveguide, 
where the field is concentrated. Despite this, we have occasionally been able to obtain 
tuning ranges of ^6 nm. 

To measure the change in lifetime when the NV center becomes coupled with the 
cavity, we first tune the cavity resonance using plasma etching until it is to the blue 
of the NV center’s ZPL. We then gradually red tune the cavity onto resonance with 
the ZPL while measuring the lifetime using a filtered super continuum laser (NKT 
photonics, SuperK EXTREME) as a pulsed laser source around 532 nm. As shown 
in Eigure |1.15[ the lifetime decreases from 18.5 ns when the cavity is completely off 
resonant to 11.6 ns when the cavity is tuned to the ZPL. This gives a Purcell factor of 
P = 0.59, or a ZPL enhancement of P' = 20 for ^ = 0.03. We emphasize here that this 
enhancement factor implies that the ZPL emission is now ^40% of the total emission 
from the NV center. 


1.10 Outlook 

Eollowing the experiments presented in Sectiou pTBl there has been tremendous progress 
in using optical photons to connect remote NV centers. Euture work will explore the 
possibility of combining the few-qubit local spin register with the resource of remote 
entanglement through, for example, teleportation of quantum information between 
remote, long-lived nuclear spin memories. Distribution of quantum information by 
teleportation has been demonstrated in several different physical systems (Nolleke 


et a/., 2013 Steffen et a/., 2013 Krauter et a/., 2013). However, remote teleportation 


over long distances is an outstanding challenge with solid state systems, which, in the 
long run, offers the possibility of implementing more scalable local quantum nodes. 
In particular, one major factor in determining the ultimate bandwidth and fidelity of 
a remote quantum channel is the efficiency with which single photons can be made 
to interact with the stationary qubits and then extracted into, for example, a optical 
fiber. Current photon collection efficiencies and entanglement rates may be sufficient 
for proof-of-principle experiments, but certainly must be greatly improved for the im¬ 


plementation of practical quantum networks. The work described in Sections |1.8| and 
1.9| begins to address this issue. The spectral stability of NV centers is an important 
factor determining the repetition rate of experiments requiring indistinguishable op- 
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Fig. 1.15 Purcell enhancement of the NV center ZPL emission using hybrid photonic crystal 
cavities, a. SEM image of hybrid diamond-PMMA cavities. Visible are the top surface of 
the suspended triangular beams, the top of the the PMMA layer and slabs, along with the 
exposed diamond surface underneath, b. Transmission spectrum of hybrid cavity (magenta) 
and emission spectrum of NV center (black). The spectra have been scaled relative to each 
other to fit on the same scale, c. Zoomed in spectra showing ZPL of NV center (black) and 
cavity resonance tuned to three different frequencies (blue, magenta, and yellow), d. Lifetime 
measurements corresponding to the detunings in c., taken by exciting the NV center with a 
pulsed laser and collecting the PSB emission. See (de Leon et al, 2014) 


tical photons, such as remote entanglement and teleportation. The availability of NV 
centers with lifetime-limited linewidths would eliminate the need for active stabiliza¬ 
tion and post-selection of the the transition frequency in future experiments. Cavity 
QED has been used to enhance the interaction of photons with atomic systems and 
improve their collection efficiency, leading to demonstrations of quantum teleportation 
with high fidelity and success rates (Nolleke et a/., 2013). Diamond-based nanopho¬ 
tonics offers the possibility of achieving efficient photonic quantum devices in a solid 
state system. 

Finally, future quantum technologies will most likely be implemented with a combi¬ 
nation of several types of physical systems, each with its own strengths and weaknesses. 
There are already many theoretical proposals and experimental demonstrations of in¬ 
terfacing NV centers with other quantum devices (Kubo et a/., 2011 Rabl et a/., 2010). 
Most of these demonstrations make use of the NV center’s unique spin properties, ef¬ 
fectively turning its interaction with the magnetic environment from a detriment into 
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a resource. Analogous to this, an potential direction for future exploration would be to 
take advantage of the NV center’s optical transitions to interface it with other systems. 
Combined with the ability to communicate information with remote systems through 
optical photons, the NV center becomes a quantum system that not only has uniquely 
attractive properties as an isolated system, but can also be made to interact with other 
systems in a controlled manner. Therefore, the ideas and techniques described in this 
chapter will likely enable researchers to use the optical properties of the NV center as 
a powerful resource for building complex and useful quantum technologies. 
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